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Bucher jlj has recently proposed an interesting brane-world cosmological scenario where the “Big 
Bang” hypersurface is the locus of collision of two vacuum bubbles which nucleate in a five di¬ 
mensional flat space. This gives rise to an open universe, where the curvature can be very small 
provided that d/Ro is sufficiently large. Here, d is the distance between bubbles and Ro is their size 
at the time of nucleation. Quantum fluctuations develop on the bubbles as they expand towards 
each other, and these in turn imprint cosmological perturbations on the initial hypersurface. We 
present a simple formalism for calculating the spectrum of such perturbations and their subsequent 
evolution. We conclude that, unfortunately, the spectrum is very tilted, with spectral index n s = 3. 
The amplitude of fluctuations at horizon crossing is given by < (Sp/p) 2 >~ (Ro/d) 2 Sfdk 2 , where 
Se ^ 1 is the Euclidean action of the instanton describing the nucleation of a bubble and k is the 
wavenumber in units of the curvature scale. The spectrum peaks on the smallest possible relevant 
scale, whose wave-number is given by k ~ d/Ro. We comment on the possible extension of our for¬ 
malism to more general situations where a Big Bang is ignited through the collision of 4D extended 
objects. 

I. INTRODUCTION 

The possibility of having gravity localized j||| on a 
brane which moves in a higher dimensional space has re¬ 
cently stimulated the search for alternatives to the stan¬ 
dard inflationary paradigm |^,|l|,||. In these alternative 
scenarios, the Big Bang would result from the collision 
of 4-dimensional extended objects propagating in five di¬ 
mensions. Pioneering work in this direction Q did not 
use the localization of gravity as an essential ingredient, 
but this was at the expense of a somewhat singular be¬ 
haviour at the moment of collision when the “bulk” 
in which the branes propagate momentarily disappears. 

The generic predictions of inflation seem to be in good 
agreement with current cosmological data, and therefore 
any alternative proposal has to measure up to high stan¬ 
dards. In particular, it is important to clarify the mecha¬ 
nisms by which cosmological perturbations are generated, 
since at the moment these provide the finest tests for any 
theory of initial conditions. The purpose of this paper is 
to consider this problem in the context of the model pro¬ 
posed in Ref. [jlj, where it seems to be best posed (see 
also §). 

The paper is organized as follows. In Section II we 
briefly review Bucher’s scenario, describing the geometry 
of the problem. In Section III we derive the expression 
for the primordial cosmological perturbations in terms of 
the physical parameters of the model. In Section IV we 
discuss the subsequent evolution of such perturbations, 
and in Section V we summarize our conclusions. 


v = dR/dT = tanh (3. 


II. A BIG BANG FROM BUBBLE COLLISION 

In Bucher’s model one starts with a metastable (or 
“false”) vacuum in five dimensions, which is flat space 
or a very mildly expanding de Sitter space. This decays 
through bubble nucleation into an anti-de Sitter (AdS) 
phase, where the five dimensional cosmological constant 
is negative. The model requires the existence of degen¬ 
erate discrete AdS vacua. Upon collision of two bubbles 
corresponding to different vacua, a domain wall forms, 
where gravity is localized a la Randall-Sundrum ||. The 
“Big Bang” hypersurface is the locus where the world- 
sheets of the two bubbles meet, and the domain wall, or 
“local brane”, is the place where we are supposed to live. 

Let us use coordinates X A = ( X l , W, T) in the original 
5D Minkowski space, where a pair of bubbles of radius 
R 0 nucleate at T = 0, separated by a distance 2d. The 
radius Rq is a fixed parameter of the theory, which is 
related to the bubble wall tension a and to the energy 
gap e between the original Minkowski phase and the final 
AdS phases || (this relation, however, will not be needed 
in the following discussion.) Let us choose the origin of 
coordinates to be at the center of one of the bubbles. 
After nucleation, the bubble wall expands with constant 
acceleration, following a hyperbolic trajectory which can 
be parametrized as 

T = Ro sinh (3, (1) 

R = R 0 cosh (3, (2) 

where R = (|A| 2 + IV 2 ) 1 / 2 , and (3 is the boost parameter 
of the bubble wall 
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The collision takes place on the plane W = d, along the 
hyperboloid 

T 2 - |Af = d 2 - R 2 0 . (3) 

To the future of this hyperboloid, and on the plane 
W = d, a domain wall forms where gravity is four di¬ 
mensional [|||]]. Throughout this paper, we shall asurne 
that, locally, the collision process is almost instantaneous 
compared with the lengthscales of our interest, and that 
a fixed fraction of the energy of the collision is channeled 
into the degrees of freedom which live on the local brane. 
We shall also assume the usual fine tuning between the 
tension of the local brane and the AdS radius, so that the 
effective 4D cosmological constant is sufficiently small. 

In the unperturbed setup, the 4D Minkowski space on 
the plane W = d is matched into an open FRW model 
along the hypersurface (|j). Along this hyperboloid, the 
metric is continuous, although its derivative is not (this 
jump in the extrinsic curvature corresponds to the jump 
in the Hubble rate accross the surface). It is convenient 
to use Milne coordinates to describe the 4D Minkowski 
space before the collision, 

^ds 2 = — dt 2 + t 2 dXl 2 H , (4) 

where the t = const, surfaces have the hyperbolic geom¬ 
etry described by d£l 2 H . In terms of Minkowski coordi¬ 
nates, Milne time is given by t = ( T 2 — \X\ 2 ) 1 ^ 2 . This 
form of the metric is valid for t < Rq cosh -1 (d/iio)- 

For t > Rq cosh _1 (ci/I?o), the 4D metric on the “local 
brane” is given by 

^ds 2 = — dt 2 + a 2 (t)dVl 2 H . 

The local brane separates the interiors of the two vac¬ 
uum bubbles, which have already collided and continue 
to expand. Due to the 0(3,1) symmetry of the pro¬ 
cess, the bulk metric on both sides of the local brane is 
given by Schwarzschild-AdS § 0 - The evolution of the 
scale factor after collision is given by the usual Friedmann 
equation 


Here, H = a/ a, n 2 = 8nG is the effective 4D gravita¬ 
tional coupling and p is the matter energy density which 
is deposited on the local brane after collision. This will 
also include some “dark radiation” Q, because a signif¬ 
icant fraction of the energy of the colliding bubbles may 
not stick to the brane; it may just fly into the bulk con¬ 
tributing to the effective mass of the Schwarzschild-AdS. 
In (||) we have also neglected “brany” corrections in the 
right hand side, proportional to K 4 p 2 £ 2 . These will be 
unimportant provided that the AdS radius t in the bulk 
is sufficiently small. The discussion of cosmological per¬ 
turbations with these corrections included is significantly 
more complicated and it is left for future research. 


From the spatial curvature scale of the open uni¬ 
verse at the moment of collision is given by (d 2 — Rq) 1 / 2 , 
so the flatness problem is solved provided that d is suf¬ 
ficiently large (as discussed in Section III.A, this gives a 
lower bound on d/Ro). Of course, the homogeneity prob¬ 
lem is also solved, due to the residual 0(3,1) symmetry 
of the colliding bubble setup along the hyperboloid (;3). 
Thus, in principle, Bucher’s scenario seems to provide an 
attractive starting point for a brane-world cosmology. 

III. INITIAL PERTURBATIONS 

Aside from solving the homogeneity and flatness prob¬ 
lems, a serious candidate for a theory of initial conditions 
should explain the origin of the primordial cosmological 
perturbations. In brane collision scenarios these may be 
seeded by preexisting fluctuations in the shape of the 
colliding branes. These produce distortions of the “Big 
Bang” hypersurface as well as initial perturbations in the 
distribution of the energy density. In this Section we es¬ 
timate this effect for the case of colliding bubbles. 

Quantum fluctuations on expanding vacuum bubbles 
have been studied in Refs. |ll[. For a bubble which 
is centered at the origin of coordinates, the perturbed 
worldsheet X A is conveniently parametrized as 

x A (e) = x A (n+N A (nx(n = a+ x /r 0 )x a . 

Here, N A is the unit normal to the unperturbed world- 
sheet of the bubble, which has the internal geometry of a 
4D de Sitter space of radius Rq , are a set of coordinates 
in this space and X A stands for the unperturbed world- 
sheet (Q). It is known Jll]] that the normal displacement 
X behaves like a worldsheet scalar field with the tachy- 
onic mass m 2 = —4R^ 2 , which obeys an equation of the 
form f] 

- □* - 4Rq 2 x = 0 . ( 6 ) 

Here, □ is the covariant d’Alembertian in a 4D de Sitter 
space. 

We must also consider the second bubble, whose cen¬ 
ter is at a distance 2d from the first. For perturbations 
which are symmetric with respect to the plane W = d, 
the collision will still take place on this plane, the sec¬ 
ond bubble being just a mirror image of the first. In 


*Here, and in the following discussion, we are neglecting the 
self-gravity effect of the bubble walls on the perturbations. 
Domain walls produce a repulsive constant gravitational force 
which causes an acceleration of order (Gct) - 1 , where G is the 
five dimensional Newton’s constant and a is the wall tension. 
If this length is very large compared with the inverse of the 
proper acceleration of the bubble wall, given by Rq 1 , then 
the self gravity of the wall will be negligible. 
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what follows, we shall restrict our attention to such “Z 2 - 
symmetric”perturbations. It is easy to show, following 
a calculation similar to the one presented below, that 
the antisymmetric mode has an amplitude comparable to 
that of the symmetric mode, but does not contribute to 
cosmological perturbations to linear order, so its effect 
seems to be much smaller than that of the symmetric 
mode. 

An important observation is that, by conservation of 
momentum, the fluid lines will be orthogonal to the 
perturbed surface of collision, which is therefore a co¬ 
moving, or velocity orthogonal surface 1121. It is known 
(see next section) that the growing mode of the curvature 
perturbation 1Z C on co-moving surfaces is a constant of 
motion on scales larger than the Hubble rate, and hence 
1Z C will be a quantity of interest to us. In the present 
case, however, it will be equally important to consider 
the decaying mode, since its initial amplitude is compa¬ 
rable to that of the growing mode. Therefore in order 
to characterize the initial perturbation we will also need 
the density perturbation S c on the same initial surface. 


A. Initial value of 1Z C 

Due to local shifts St in the time of collision caused 
by fluctuations of the bubble shape, the Big Bang sur¬ 
face will no longer be the smooth hyperboloid (|). As 
explained in the preceeding paragraphs, we can restrict 
attention to perturbations which are symmetric with re¬ 
spect to the plane W = d. Then, to linear order in St, the 
metric induced on the Big Bang surface is easily obtained 
from (^), and it is given by 

^ds 2 = (l + 2 S -pjt 2 dtf H . 

Since we live in an almost flat universe, we shall be in¬ 
terested in lengthscales which are short compared to the 
curvature scale of the unperturbed surface. On these 
scales dfl H can be replaced with a flat metric. Then, the 
intrinsic curvature R of the perturbed surface is easily 
found to be 


{3) R= —ATI - -A— 
" t 2 c ~t 2 t' 


where A is the co-moving Laplacian. The first equality 
is just the conventional definition of the curvature per¬ 
turbation used by most authors Jl6| . In a flat FRW with 
arbitrary scale factor and in an arbitrary gauge this def¬ 
inition takes the form 


i3) R = ^AR, 
a z 


(7) 


To proceed, we must find St as a function of the normal 
displacement y. Since all points on the unperturbed hy¬ 
perboloid are equivalent, we may just consider the vicin¬ 
ity of the point X = 0, where t = T. It should be noted 
that the field \ lives on the unperturbed bubble, which 
hits the collision plane at some value Pd of the unper¬ 
turbed boost parameter determined by d = Rq cosh pd- 
However, the actual collision does not take place at that 
time, but at the time when the perturbed bubble reaches 
W = d, at some value (3 C of the unperturbed boost pa¬ 
rameter given by d = {Rq + x) cosh (3 C . Thus, the shift in 
the time of collision is given by 


St = ST = T-T 

= (R 0 + x) sinh f3 c - R 0 sinh/3 d « — . } . 

smh / 3 d 


We conclude that the curvature perturbation on the ini¬ 
tial co-moving surface is given by 


_ ST 
R c {i) = — = — 


X 


T 


Ro sinlr I3 d 


( 8 ) 


Let us now estimate the size of this effect. 

As mentioned above, x satisfies the equation for a 
scalar field on the worldsheet de Sitter space, Eq. §)■ 
The corresponding canonical field p with dimensions of 
mass is related to x by fll|] 

</> = o- 1/2 x, 


where a is the tension of the wall. As is well known, when 
a given mode of a nearly massless scalar field crosses the 
de Sitter horizon Rq 1 , it “freezes” with some amplitude 
4>k ~ Rq 1 ■ However, in our case this amplitude does not 
stay constant. Due to the tachyonic mass m 2 = —AH 2 
in (0), a mode with wavenumber k grows proportionally 
to er as 


Xk{P) ~a- 1/2 I? 0 -V-^, (9) 

where f3k is the value of the time-like boost parameter 
[see Eq. (|J)] at which the physical wavelength of the 
mode k~ 1 R( j e >3 becomes larger than the inverse expansion 
rate on the bubble, Rq. Here we adopt the convention 
that the co-moving curvature scale corresponds to k ~ 1. 
It is easy to see that at the time of collision we have 
exp(/3 d — Pk) ~ k~ l {d/ Rq), and therefore, 

Xk(Pc) ~ (J ^ 1,2 -^ k ~ 1 - 

Using exp Pd ~ d/Ro in Eq. (|j) we find that by order of 
magnitude, 


f* -1 - (10) 

Here Se ~ <jRq 1 is the Euclidean action of the in- 
stanton describing the nucleation of the bubble. 


where now ^ R is the perturbation of the intrinsic cur¬ 
vature scalar in the constant time surfaces of the corre¬ 
sponding gauge. 
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This perturbation is rather minute even on scales com¬ 
parable to the curvature scale k ~ 1. An upper bound on 
Ro/d can be obtained as follows. The spatial curvature 
radius of our present universe is much larger than the 
present Hubble radius Hq 1 . Taking into account that 
the curvature radius of the initial surface was given by d, 
this leads to the constraint 

Zid > H 0 , 


is the derivative of the perturbation with respect to the 
proper time r = R 0 /3 measured by an observer on the 
bubble. Thus, the change in the boost parameter is given 
to linear order by 

&P = P - Pd = X- TT coth Pd- 

Kq 

With these relations, we obtain the density perturbation 
at the moment of collision as 


where Z{ is the redshift at which the collision surface is. 
This redshift is given by 


Zi 




1/4 


where we assume radiation dominance for most of the 
cosmological evolution. In the last inequality we have 
used the fact that the energy density after collision will 
not be larger than the input energy density, which is of 
order a cosh (3d- From the two previous inequalities we 
find 


§ < icr 24 s 


1/5 

E 


Therefore, we have 

n c(i) (k) <io~ 24 s- 3/w k-\ (ii) 

If there was no initial density perturbation, then this 
quantity would stay nearly constant until horizon reen¬ 
try (see next section). This would lead to an amplitude 
of perturbations smaller than the observed value by some 
20 orders of magnitude. As we shall see, the bound 0) 
still applies to scales comparable to the curvature scale. 
However, on smaller scales the initial density perturba¬ 
tion <5 c (j) induces a much larger perturbation on 1Z C by 
the time of horizon reentry. Unfortunately, this comes at 
the prize of a very tilted spectrum, which is not compat¬ 
ible with observations. 


B. Initial value of 8 C 

To evolve the cosmological perturbations, two initial 
conditions are needed for each wavelength. In addition 
to the initial curvature 1Z C we will also need the density 
perturbation S c . By assumption, a fixed fraction of the 
energy of the collision goes into the brane, and therefore 
the density at the moment of collision is proportional to 
the Lorentz factor 7 = cosh/3. Since the boost parameter 
is additive, the perturbed one will be given by 

P = Pc + X; 


where 


J_<9x 

R 0 dp 


5 c u) = tanh PdSp = xtanh j3 d - ( 12 ) 

ito 

As noted above, in the case of our interest, the r.m.s. fluc¬ 
tuation in x grows exponentially fast with /3. However, 
the combination which enters S c anihilates the leading 
term in this exponential dependence, and leaves only a 
contribution which decays with [3. From this argument 
alone, it should be clear that the spectrum of 5 C will not 
be scale invariant, and that it will have opposite tilt to 
the spectrum of 1Z C . To extract the spectral index, we 
need to look at the detailed form of the mode functions. 

For definiteness, we shall use the open chart on the de 
Sitter worldsheet of the bubble. The conclusion, however, 
is independent of our slicing since we are interested only 
in wavelenths much smaller than the curvature scale. In 
this chart, equation (||) for the evolution of x reads 

3 4 1 + fc 2 

Xk + — coth [3xk - -hiX fc + p2 . ,2 a * k = °' ( 13 ) 

-rto -fi 0 Ay sum p 


The quantum state of perturbations on a nucleating bub¬ 
ble is uniquely determined by de Sitter invariance [ |l3| . 
It is given by the so called Bunch-Davies vacuum (also 
known as Euclidean vacuum). The corresponding modes 
in the open chart have been studied by jl4 151, and for a 
scalar field with mass m 2 = — 4Pq they are given by 


Xk(P) oc 


P 2 lfe ( cosh/3) 
sinh /3 ’ 


where P 2 fe are the Legendre functions with the branch cut 
from —00 to 1 on the real axis. In the above equation we 
have ignored a contribution proportional to the decay¬ 
ing mode, which is accompanied by the factor e~ nk and 
which is therefore irrelevant at wavelenths much shorter 
than the curvature scale. Expanding the Legendre func¬ 
tion for large p, we have 


Xk(P) 


sinh /3 fe 


sinh /3 + 


k 2 + 4 
6 sinh /3 


+ O(e~ 30 ) 


(14) 


Here, the normalization is fixed as in Eq.(^), the only 
difference being that we have been careful to keep more 
terms in the expansion because the leading one cle arly 
does not contribute to ©• Substituting 0 ) in ( flq ) 
and using exp /3 ss d/Ro we have 


dc(i) 
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(15) 










For k ~ 1, this is of the same order as the curvature 
perturbation given in where both are very small. 

However, © can be much larger on small scales. Un¬ 
fortunately, this is due to a strong tilt in the power spec¬ 
trum, corresponding to n s = 3 in the standard notation. 


IV. EVOLUTION OF THE PERTURBATIONS 


Eqs. (j 8 j) and © provide the initial conditions (po¬ 
sition and momentum, as it were) for the evolution of 
cosmological perturbations in the FRW phase. As men¬ 
tioned above, since our universe is reasonably flat, the 
scales of interest to us will be much smaller than the cur¬ 
vature scale and it is a good approximation to consider 
perturbations to a flat FRW universe. 

In an arbitrary gauge, the perturbed metric for scalar 
perturbations can be written as |l(| 


ds 2 = a 2 (rj) 2 


(1 + 2 AY)drj 2 - 2BY i dgdx i 


+ [(1 + 2 DY)S i:i + 2EYij] dx'dx 0 


where, Y oc e* kx is the appropriately normalized plane 
wave, and summation over modes with different wave 
number k is omitted. The vector and the traceless ten¬ 
sor constructed from Y are defined by Y, = — and 

Y/j = k~ 2 Y.ij + \8ijY. We also use convenient combina¬ 
tions of metric perturbations defined by 

1Z = D + -E, 

ka g = E' — kB , (16) 

where a prime denotes differentiation with respect to rj. 
1Z is the quantity related to the perturbation of the spa¬ 
tial scalar curvature of a time slice through Eq. ( 0 ), 
and <T g is related to the shear of the hypersurface normal 
vector field | The familiar gravitational potential in the 
Newton gauge is gauge invariantly defined as 

& = -K + Hk- 1 <r g , (17) 

where 7 i = a'/a. 

The energy momentum tensor of a perturbed perfect 
fluid is given by 

= (p + P)vf l u v + Pg^, (18) 

with p = p + YSp , P = P + Y5P, u° = a _1 ( 1 — AY), 
and u 1 = a~ 1 vY l . We shall assume that the fluid con¬ 
sists of a single component. Then the ratio between the 


Wor the present discussion, the geometric interpretation of 
a g will be irrelevant, it will just be used as a convenient 
variable. 


density and pressure perturbations becomes a function of 
the background energy density. We denote this ratio by 
c 2 = SP/Sp. 

In the comoving gauge, in which v — B = 0, the per¬ 
turbed Einstein equations become 


HA = 7 Z’ c , 

(19) 

2 2 

di(kag) k 2 lZ c = 2 Sp, 

( 20 ) 

(ka g y + 2 H(ka g ) - k 2 + iz)j = 0 , 

( 21 ) 

w '“” 2 

( 22 ) 


Eqs.©), ©) and (H^) have already been simplified by 
using (|ig|). For reference, we also quote the background 
equations 

K 2 n 2 K 2 n 2 

H' = - — {p + ZP), n 2 = —p. (23) 

Eqs. ©) and (||^) give the following expression for the 
gauge invariant $ in terms of the density perturbation 5 C 

The value of $ is larger than 8 C by the factor of 7d 2 //c 2 , 
and hence this is not a particularly illuminating variable 
on very large scales. 

The evolution equations for given initial values for 1Z C 
and S c can be obtained by combining Eqs(|o]), (j|I|) and 


dlZ c 2 5c 
In = Cs (1 + re) ’ 

^-^-(a 2 7fd c ) = fc 2 ( 1 + w)7Z c , 
a z dN 


(25) 

(26) 


where w = P/p is the parameter characterizing the equa¬ 
tion of state and dN = Tidrj. 

To discuss the evolution of perturbations it is useful to 
formally integrate the second equation, 


1 r N k 2 

5 ^^HjJ“ 2 ' H) W ll+w)K ' dN 

*>• (27) 

The function a 2 /Ti is an increasing function, so pro¬ 
vided that 1Z C stays approximately constant at late times 
(which will be a self-consistent assumption), the inte¬ 
gral in the right hand side is dominated by the contri¬ 
bution from the neighborhood of the upper boundary 
of integration. If we consider a simple case in which 
c 2 = w =const., the scale factor and H are given by 
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« = (v/vo ) 1+3 " , U = , . o V 

1 + 3w 


Performing the integral in 
- 2 + 2 w k 2 


we have 


5 + 3 w K 2 


TZ C + 


a 2 Ti 


(fl 2 Wc) (i) , 


(28) 


(29) 


where we have assumed 1Z C ~ const, at late times [i] Ini¬ 
tially, 7 Z c 1$ 6 C and consequently the first term in the 
right hand side is suppressed with respect to the second 
at least by a factor of (k?/TL 2 ) <C 1. Substituting the 
dominant part into Eq.(|25|), we have 

Kc(f) ~ K*) ~ (a ktS)^ j \j rW - (30) 

Then, this integration is performed to obtain 

'R-c(f) ~ ~ 3(1 _ u; 2 ^ c (*)‘ (31) 

In obtaining m we have neglected the contribution from 
the first term in the right hand side of (|29|). It is easy 
to check that this term does not have any effect until 
the wavelength of the mode is comparable to the bubble 
radius k ~ Tt. At that time, the second term in (^) is 
unimportant, and the amplitude of density perturbations 
at horizon crossing can be read off from the first term 


2 + 2w k ^ 

5 + 3 wTP c(/) ' 


(32) 


Finally, the Newtonian potential at horizon crossing can 
be found from the relation (pi|) 


Substituting ( 31]) into the previous two equations and 
using © andp^) we have 


*\hc ^ ^c|/ic 


n c(f) 


Ro 0 - 1/2 
d bE 


k, 


which of course displays the same strong spectral tilt as 
the initial density perturbation S c ^y 


V. CONCLUSIONS 

In this paper, we have discussed the primordial spec¬ 
trum of density perturbations in the brane-world model 
proposed by Bucher JTJ . In this scenario, two bubbles 


^Incidentally, even when TZ C is not constant, but increa sing 
with time, the order of magnitude of the first term in ( |2S| ) 
does not change. 


nucleate in the 5-dimensional Minkowski bulk, and their 
collision forms a brane where gravity is localized. This 
model solves the homogeneity and flatness problems, pro¬ 
vided that the separation 2d of the bubble nucleation 
points is sufficiently large compared with the bubble ra¬ 
dius Rq. 

We have evaluated the initial spectrum of scalar den¬ 
sity perturbations, assuming that they originate from 
quantum fluctuations on the bubbles. An important ob¬ 
servation is that the surface of collision coincides with 
the so-called comoving hypersurface (or velocity orthog¬ 
onal slicing). On this initial surface, we first evaluated 
the density contrast 6 C = Sp/p and the spatial curvature 
perturbation 1Z C [See Eq. ( |l6| ) for definition]. The ini¬ 
tial value <5 c (j) turns out to have a very steep spectrum 
corresponding to n s = 3, while the spectral index for 
lZ r (i) is n s = — 1. The amplitude of both perturbations 
is comparable at wavelengths of the order of the curva¬ 
ture scale, and therefore 6 C ^ has a larger amplitude on 
scales relevant to present observations. 

It is known that the curvature perturbation in the co¬ 
moving gauge, 1Z C , is conserved on scales much larger 
than the horizon, and therefore this quantity is often used 
to discuss the evolution of perturbations in the early uni¬ 
verse. The constancy of 1Z C , however, does not hold for 
the decaying mode. In the present case, the initial con¬ 
ditions which arise as a result of bubble collision, contain 
a significant amount of “contamination” from the decay¬ 
ing mode. Therefore, 1Z C does not stay constant in the 
subsequent evolution of perturbations, which is described 
by two coupled first order differential equations. Given 
initial values 6 C ^ and TZ c (i ), we solved these equations 
to find that the final value of 1Z C (at the time of hori¬ 
zon crossing) becomes comparable to the initial value of 
the density perturbation 5 C ^ [See Eq. (|3^)]. This final 
value of TZ C also gives the order of magnitude of the den¬ 
sity contrast at horizon crossing. Thus, we conclude that 
the spectral index for primordial density perturbations is 
n s = 3. The estimated amplitude is very small at wave¬ 
lengths comparable to the curvature scale. Perturbations 
of 0(1CT 5 ) at the present horizon scale may of course 
be obtained by choosing the ratio between the curvature 
scale and the present horizon scale appropriately, but the 
spectrum is too steep to be consistent with observations. 

Although in this paper we have investigated a par¬ 
ticular realization of the brane Big Bang, we may have 
learned a few lessons which may be useful in more general 
cases. First, in the present model the bubble fluctuations 
are described by an effective 4-dimensional scalar field 
with a negative mass squared. Consequently, fluctuations 
on the bubble worldsheet have a red spectrum, which has 
a larger amplitude for longer wavelengths. However, the 
resulting spectrum of density perturbations turned out 
to be a blue one. This means that in principle it may be 
possible (although perhaps not easy) to generate a nearly 
scale invariant spectrum even when the mass of the ef¬ 
fective field corresponding to the bubble fluctuations is 
not close to zero. Second, we need to be careful in using 
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“standard” results of cosmological perturbation theory, 
which in some cases neglect the contribution from the 
decaying mode. In order to determine the evolution just 
after collision, two initial conditions must be supplied, 
and both turn out to be important. 
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